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Technical Considerations and Recommendations for SELENE 


Originally conceived for the purpose of delivering power to a lunar station during the 14-day 
lunar night, the SELENE program during the past year evolved into a far more significant one for long- 
term NASA interests as it became apparent that the system could also 

• supply the power for the orbital transfer vehicles needed to economically transfer materials 
and personnel from LEO to space platforms, lunar stations, and beyond, both to establish them 
in the first place and thereafter to sustain them, and 

• augment or replace new satellite power and propulsion systems and regenerate aging high- 
value satellites. 



Fig. 1 Loser Power Beaming Elements 


The primary elements of a laser power beaming system required to satisfy these applications — surface, 
transportation., and satellite - are shown in Fig. 1. Although the elements can be separately identified 
and developed, there is a strong interplay between them, and this report will attempt to look at them as a 
whole and to make recommendations for the future. The views expressed here are those of the author, 
and do not necessarily represent those of NASA or of any of its other contractors. 


The two most important factors 
driving the configuration of the power 
beaming system have been (1) minimization 
of weight in space because of the extremely 
high costs associated with it and (2) the 
continued use of photovoltaic conversion with 
solar cells so that the receivers can function 
with both solar and laser power. The spectral 
response characteristics for the major types 
of available cells are shown in Fig. 2 from 
which it can be seen that the high efficiencies 
required to minimize weight in space will 
recommend operation with GaAs in the 0.8 - 
0.85 p range. Silicon is the most widely-used 
material today, and it is expected that 
radiation damage will move its response peak 
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into that region as well. Technical difficulties with the laser and optics become dramatically worse at 
shorter wavelengths, so AIGaAs offers no advantages, and the low efficiencies of the IR cells would 
offset their other benefits (although ocular damage could ultimately prove to be an overriding 
consideration if beaming were to take place below the Earth's limb.) 

Low power semiconductor and Ti:sapphire solid-state lasers are the only conventional devices to 
meet this unique wavelength requirement today, but the free electron laser (FEL) at least in concept 
offers virtually unlimited wavelength flexibility together with the potential of reliable, extremely high power 
operation without the massive flows associated with fluid mediums or the optical degradation of solid 
ones. Because of its demonstrated high power capability (over 1 GW at microwave frequencies) and the 
potential reliability of its all-solid-state and highly-modular architecture, the induction linear accelerator 
FEL was the system of choice during much of the SELENE program. The induction linac uses 
independent cavities which each accelerate the electron beam through a fixed voltage via magnetic 
induction and, being a low impedance machine, can accelerate very high currents in pulses with lengths 
on the order of tens of nanoseconds at several kilohertz repetition rates. However, these characteristics 
that made the induction linac so attractive in the past as the beam source for a weapon have proved to 
be its undoing for laser power beaming to photovoltaic receivers, and interest has turned to the use of 
other types of accelerators, specifically the RF linac and the Microtron. The RF accelerator, which unlike 
the induction has already been operated at short wavelengths, uses a microwave field with a highly- 
relativistic phase velocity to accelerate the electrons, and, though limited to much lower peak currents, 
produces closely-spaced micropulses with lengths of tens of picoseconds, thereby producing quasi-CW 
interactions. From the operational standpoint, they and their klystron power supplies are proving to be 
very reliable, and years of valuable experience have been gained with them at SLAC as well as at other 
facilities around the world. Another very promising concept, the Microtron in which the electron beam 
travels in a tightly-spiraled path in a fixed magnetic field and is recirculated through a set of RF cavities, 
has been extensively pursued by investigators in the former Soviet Union, and collaboration on it is now 
underway in the US. 

The free electron laser is the only single means of producing coherent radiation over 
frequencies extending from the microwave into the soft X-ray, and is therefore likely to continue 
to receive attention from at least the university sector. However, the types of applications of 
interest to the researchers, together with the cost of components and facilities for high power 
operation, will limit the scope of those efforts, and the R&D that will be needed for the class of 
FELs that would be suitable for power beaming, waste treatment, bulk sterilization, and the like 
will cease if funding from DoD and NASA is curtailed. The excesses of the past have arisen 
from weak management and cost controls, not from fundamental and unavoidable technical 
difficulties, and a well-run unified program needs to be defined and undertaken so that the 
nation can finally benefit from this very promising technology. 

The problems with pulsed operation of a photo- 
voltaic array can be seen in Fig. 3 which essentially 
gives the variation of conversion efficiency with total 
series resistance (load plus internal) and short-circuit 
current, the latter being proportional to a suitably aver- 
aged illumination level. While the maximum efficiency 
shows a fairly substantial logarithmic increase as the 
intensity rises, excursions of more than about a factor 
of two from optimum resistance can dramatically 
degrade performance. To put matters in perspective, 
one solar unit corresponds to about 1,^ = 210 4 with 
silicon, while the power beaming program is interested 
in pulsed irradiation with time-averaged levels of at 
least a few suns and with peak levels about 1000 times 
the average. The optimum load varies approximately 
inversely with the illumination level, and very high 
intensities can eventually cause the internal cell 
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resistance alone to exceed the optimum value. Landis suggested two years ago that high currents 
resulting from strong irradiation could produce unacceptable losses in the internal series resistance, but 
that the minority carrier lifetime of the cells (a measure of the time for recombination after light-induced 
charge separation) could act to smear shorter pulses, effectively 
integrating them over that scale and reducing the peak currents. 

The result would be such that the large values for indirect bandgap 
materials like silicon would integrate the pulses in both cases, 
causing the cells to respond to the average intensity, whereas the 
shorter lifetime for direct bandgap materials like GaAs (comparable 
to the induction pulses) would smear only the short pulses (Fig. 4) 
and cause the cells to respond to the peak intensity of the longer 
ones with high currents and large series-resistance losses. This 
was confirmed by measurements in Feb 92 at LLNL and Rocket- 
dyne by NASA personnel using conventional lasers with simulated 
induction and RF formats: the results showed that silicon operated 
nicely with either, but that both conventional and concentrator GaAs 
cells gave far better performance with the shorter pulses. 



induction ra. 


Hg. 4 GaAs Photovoltaic 
Response to FEL Irradiation 


Various means for alleviating the problem for the induction FEL were studied by this author, 
including both capacitive loading and series connection to reduce the resistance, and pulse lengthening 
via splitting and delay to reduce the peak currents, but none proved to be particularly cost-effective, and 
efforts were terminated when it was realized that the problem lay not just in the pulse duration but in the 
interpulse time as well. The fact is that the required accuracy and dynamic range may be quite difficult 
and costly to achieve if the system is to be truly multi-purpose by being be able to accommodate both 
solar irradiation and the large peak currents from laser pulses. Fortunately, the minority carrier lifetimes 
of both GaAs and silicon are such that the pulse train from an RF FEL (at least within a macropulse) will 
be effectively averaged into a CW-like output from the cells — but the interpulse times for the induction 
FEL are long even compared to the silicon lifetime so the cell output will remain pulsed with that laser 

The problem with any laser beam is further complicated by the nearly-gaussian distribution of 
intensity about its centroid at the focus and by the centroid wander produced by both system jitter and 
atmospheric turbulence. Unless the receiver were flood-loaded with a much larger beam, the spatially 
and temporally varying illumination would require active load compensation for fully optimum operation, 
and that would considerably complicate the array design. The effects of various loading schemes have 

been analyzed by considering independent receiver 
elements each with its own load resistance, a model 
which should give a reasonable representation of the 
behavior of parallel-connected cells (see Attachment G 
for an analysis of photovoltaic operation with laser 
illumination.) The maximum total power output with 
identical elemental circuits and gaussian illumination is 
achieved if the load is optimized for about 55% of the 
peak short circuit current, and Fig. 5 compares the 
distributions of output power density and included power 
for optimum fixed resistance with those for locally- 
optimized ("active") loading for a typical case with 1,^ = 
10®. The value of the fixed resistance is higher than 
optimum in the high-intensity central region and lower in 
the wings, so the efficiency is somewhat less than would 
be achieved with the locally-optimized loading (24% less 
for the sample case.) The rather nominal performance 
penalty may prove to be a very small price to pay, 
however, when one considers that the locally-optimized 
loading would have to be "active" at relatively high 
frequencies to account for centroid motion, whereas, 
although the constant-resistance "passive" load might 
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have to be variable as well to account for changes of peak intensity (particularly for space vehicles), its 
variations would be applied on a global scale at much lower frequencies. Furthermore, while 86 5% of 
the incident power is included within the 1/e 2 radius r = 2a, that same fraction of the output power for the 
case of Fig. 5 is included within r/a = 1.89 with local optimization and only 1.55 with the optimum 
constant load, thereby providing the possibility of a smaller receiver. 

Photovoltaic cells are intended for relatively uniform flux and are not well-suited to pulsed laser 
irradiation, but their performance is understood and they are well-developed for space 
operation. The temptation to seek alternative conversion schemes should be overcome, and 
means should be sought to best adapt them to laser irradiation, but it must be established early 
on that this can in fact be done without compromising overall system performance. 

A second serious consequence of the pulse format of the induction FEL is the fact that the 
physical length of a pulse is on the order of ten meters while the pulse spacing is tens of kilometers. 
Unless an extremely long path is used as a delay line to interfere successive pulses, the induction FEL 
cannot therefore be used as a resonator, and a relatively powerful master oscillator must be used at the 
appropriate wavelength. The problem was considered difficult but at least solvable with operation at 
1.06 microns with the highly-developed NdYAG or Nd:glass lasers (LLNL wanted 30M$ to develop one 
even after all the experience that had been gained with the SHIVA and NOVA ICF systems), but, as 
noted earlier, only semiconductor and Ti:sapphire lasers operate in the 0.8 - 0.85 micron region required 
by the photovoltaic receivers, and those systems have been demonstrated only at very low power levels. 
RF FEL micropulses, on the other hand are much shorter (generally less than a centimeter), and the 
spacing between them is only a few meters, so they can easily be operated as resonators without a 
separate driver. 


Progress with semiconductor and solid state lasers should be monitored for their application to 
power beaming, not just as drivers for FEL amplifiers, but rather for use as the primary laser. 
Rapid advances being made in diode pumping promise major improvements in both efficiency 
and beam quality, and, if the elusive means can be found to lock multiple cavities together in 
phase, a myriad of coupled lasers could prove most attractive. That would especially be so if 
the entire system could function as a phased array with the locking signal coming from the 
receiver ; thereby automatically compensating for disturbances along the path. 


Turning now to propagation, a coherent monochromatic beam of wavelength X from an aperture 


with diameter D is spread with a half-angle of 
about XJD by diffraction. The resulting ideal peak 
intensity and spot size are shown in Fig. 6 at 
0.84 p with a gaussian beam truncated at 1/e 2 
(the latter is near-optimum shaping for an FEL, 
but the effect is a minor one here.) Beaming to a 


realistic receiver beyond a very low orbit will 
clearly require a large aperture, and the 12 
telescope being considered for SELENE will 
produce a 5 m spot in GEO and a 50 m one on 
the lunar surface (for comparison, the deployable 
solar panels on the Hubble Space Telescope are 
each about 12 m long.) Errors in the source 
beam and index variations along the optical path 
will cause instantaneous energy spread beyond 
the ideal size, and mechanical jitter will cause 
further time-averaged spread (kID is only 70 nrad 
with a 12 m telescope at 0.84 p), and some of the 
energy may prove to be impossible to recover. 
Whereas phase aberrations over scales compar- 
able to D will affect the distribution of the central 
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lobe itself, small-scale disturbances with an RMS phase error <> and 
lateral correlation scale d « D will scatter a factor of about l-expf-^) 
of the power into a much broader beam of half-angle TJd, effectively 
producing a composite profile (Fig. 7) with an attenuated but 
otherwise undisturbed central lobe surrounded by very low-level and 
widely-spread energy. Thus, 9.1% of the power is scattered for X720 
RMS wave distortion, and this rises rapidly to 33% at W10, so it is 
imperative that errors be small. The same is true whether the system 
is for power beaming or imaging, but there is an important distinction 
between them in that, for the same degree of degradation, the central 
lobe may be sufficiently strong to provide the required resolution and 
signal strength for imaging, whereas it may be totally inefficient for 
beaming (this can prove to be particularly important in weighing the 
significance of data from adaptive optics tests.) In either case, the 
scattered energy is effectively lost. Since FELs tend to produce 
"perfect" beams (aside from correctable errors introduced in the 
external optics), the central core should be very close to the ideal size - although certainly attenuated - 
and additional allowance need only be made for centroid motion from mechanical jitter and uncorrected 
atmospheric tilt. It should be noted that thermal blooming, a major source of central core distortion with 
many other high power laser systems, does not appear to be a problem here (at least in the expanded 
beam) primarily because of the low absorption at 0.84 p and the large apertures being used, and the 
adaptive optics system should be able to accommodate minor contributions from it as long as care is 
taken to ensure that there is always sufficient crosswind to sweep away the heated air. Nevertheless, it 
has been a major problem in the past, and should be reconsidered if changes are made in any of the key 
design parameters. 

The scattered portion of the beam arises primarily from the density variations produced by 
atmospheric turbulence, and the adaptive optics system required to compensate for them is the second 
key area, besides the free electron laser, for which extensive R&D is still required. This entails several 
intimately-related elements: 

• a synthetic beacon, target return, or other source that provides a reference signal, 

• sensors that measure the source signal and optical component motion, 

• algorithms that translate the data into commands, 

• actuators that drive the components, and 

• the phase correction elements themselves. 

The effects of turbulence on a beam inclined at 
an angle © from zenith are characterized by the Fried 
coherence length 

r Q = 0.185 (*2cos© / / C n 2(h) dh) 3/5 

for a uniform plane wave and a Kolmogoroff spectrum, 
and it has been demonstrated that appreciable improve- 
ment can be realized in both imaging and beaming 
applications if the disturbed wavefront is properly cor- 
rected over dimensions of r Q or less. A representative 
spatial distribution of turbulence strength, the Hufnagel- 
Valley profile is shown in Fig. 8, and displays the high 
level of turbulence normally found very near the ground, 
the “spike" associated with the tropospheric wind shear 
layer, and the otherwise pronounced decay with 
altitude; not shown is the highly-stratified behavior that 
arises from velocity shear and local disturbances, but 
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those variations are smoothed by the path integration for ground-to-space propagation. With the H-V 
profile, 73% of the contribution to r Q arises from the ground layer alone and only 9% is due to the "spike". 

There has been a tendency to concentrate on the path through the free atmosphere, but an 
area that should, perhaps, be of greater concern is the turbulence produced by the beam 
director itself and by thermal blooming in the high-intensity optical path within it. The severe 
turbulence near the ground is due to the earth's boundary layer, obstructions in the wind, and 
thermal gradients, so, since the power beam is directed upwards (rarely exceeding 65 degrees 
from zenith), it may be possible to design the telescope enclosure and the immediate 
surroundings to mitigate the generation of strong turbulence, thereby substantially increasing 
the turbulent coherence scale r Q . Similarly, the fundamental cause of thermal blooming is 
absorption by molecules and aerosols in the optical path, and a means must be found for 
maintaining an acceptable environment within the telescope. Neither of these need involve 
massive computational efforts, but they must be done, and the/ll require careful design 
engineering with attention paid to the basic mechanisms. 



Fig. 9 shows the results of a comprehensive series of star trail measurements taken at WSMR. 
These are unusual in that they have been presented - as all turbulence data should be-in terms of 
probability of occurrence, and this is especially important for laser power beaming to space because the 
worst-case situations of large angles from zenith and high-slew rates for which compensation is most 
difficult are also likely to be relatively rare and unduly stressing and should not be used alone to define 
the system requirements (e g., r 0 at © = 75° is 2.25 times smaller than that at small angles.) The data 
suggests that r 0 in all cases is greater than 30 mm 90% of the time; this translates to 5 cm when scaled 
from the mid-visible (0.55 p) to 0.84 p (the H-V model predicts 9 cm at the longer wavelength.) 

The design value of 3 cm dictated for use as the segment size in the SELENE program is an 
admirable goal, and should be sought if it can be done, but it is smaller than necessary for good 
overall system performance, and the penalties for using such a small size (e.g., 160,000 
elements with a 12 m telescope) should be carefully weighed. 

A complete methodology for estimating the actuator requirements for adaptive optics correction 
of global tilt and residual piston error arising from Kolmogoroff atmospheric turbulence is presented in 
Attachment E. The primary advantage of global tilt removal is that large control displacements are 
required only at relatively low frequencies, whereas, if not done, the high-frequency controls for the 
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individual segments would have to be capable of considerably larger strokes; the cost/benefits of a 
separate global tilt correction system would depend upon the particular implementation and should 
separately be addressed. Assuming an r c = 4 cm with segments of the same size in view of the afore- 
mentioned arguments, the analysis predicts RMS requirements for global tilt, piston, and subaperture tilt 
correction of 3.30 prad, 42.4 rad, and 8.53 jxrad respectively for a 12 m aperture. Using a 2.5X multiplier 
as justified in the attachment, these would translate into peak-to-peak values of 4. 1 p for edge stroke of a 
1 m global tilt corrector (12X magnification), 13.5 p for piston, and 0.05 p for additional edge stroke of 
the piston to correct for subaperture tilt. The required frequency response of the actuators is equally 
important, and that depends upon the speed at which the turbulence is swept across the beam; that 
issue is addressed as well in the same attachment, and generalized results are given for a uniform 
crosswind, a Buft on-like profile (Fig. 8), and slew-induced "wind." Subject to representative conditions 
described there, and considering two particular cases of interest, one with beaming only to high orbits 
where the velocity would be given by the Button profile and the other for beaming to a 300 km orbit for 
which slew would be the dominant wind source, the cutoff frequencies for global tilt correction would be 
2.9 HZ and 39 HZ respectively, while those for residual piston correction would be 130 HZ for the high- 
orbit system and 1 350 HZ for the low 

Low-orbit applications clearly impose extremely severe requirements on the AO system (and, 
of course, on the slew mechanism), so a shared system with a much smaller dedicated beam 
director - which would be effective at short ranges - might prove optimum for beaming to LEO 
with far fewer high-frequency elements. For example, a 2.2 meter telescope operating at 0.84 
M will produce a 2.3 m spot at 3000 km range, and this will require only 3025 4 cm elements as 
compared to the 90000 elements in a 12 m system - all of which would have to respond at 
high frequency! Although the global tilt control of the smaller system would have to respond to 
about 180 HZ, the required throws of all of the actuators would be about 4X smaller, and the 
larger system could be designed for the less-stressing high-orbit applications alone. 

An adaptive optics system corrects for phase errors by adjusting the optical path length by an 
amount equal and opposite to that caused by the disturbances; nothing is done to correct for the less 
important (but non-negligible) degradation caused by intensity fluctuations. The most accurate approach 
for phase control is via direct measure of the end-to-end optical path length variation across the beam 
with individual interferometry at each segment to define the piston displacements required for correction 
(at least within an ambiguity of 2%), but this is probably the most difficult and costly to implement by far. 
A Hartmann sensor to measure the local slope of the wavefront at each segment is less complicated and 
has become the technique of choice for adaptive optics, but that advantage is balanced by the needs to 
convert the slope information to wavefront shape via Poisson’s equation and to establish a reference (or 
references) from which to measure the shifts. The latter would be trivial with a conventional optical 
bench, but a large power beaming telescope structure is unlikely to provide such an accurate base, and 
additional measurements will be needed to relate either the heights or references of adjacent segments 
(for example, the lightweight millimeter-wave truss design being considered for SELENE is "precise" only 
at microwave scales and is quite flexible at optical ones.) 

Although the adaptive optics system should in theory be able to compensate for deflections of 
the telescope structure as well as for the disturbances from the atmosphere, most of the Poisson 
solutions involve some type of edge matching as the basis for commanding the actuators, and it is quite 
doubtful that any of the proposed edge measurement techniques will be sufficiently accurate to prevent 
an unacceptable error buildup across the beam if used throughout to give a "relative reference." 
Displacements everywhere must be controlled to within at least A/20 as noted earlier, and probability 
requires that relative displacements over a linear path involving N measurements from an arbitrary 
reference then each be accurate to a factor of at least 1/VN of that; N would be about V( 90000/jr) = 169 
for a system with 90000 segments and X - 0.84 p, and that could therefore require a measurement 
accuracy of about 30 angstroms! On the other hand, a very promising alternative appears to be the 
deterministic solution scheme proposed last year by AMP which gives the required piston and tilt 
displacements of each segment directly without iteration and without involving the error buildup of other 
techniques, but it does rely for its success upon the ability to establish an accurate reference surface 
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from which the piston displacements can be measured (Attachment F describes the application of the 
approach to hexagonal segments.) Although the predicted edge differences can in theory be used to 
individually adjust the segments (or each be combined with the local average tilt to give the piston 
differential across the gap), the differences are likely to be so small that measurement errors will 
accumulate excessively, and they should be used instead in groups to specify piston displacements from 
a common reference piston — not each with reference to a nearest neighbor There appears to be no 
obvious way of providing a single absolute reference over the scale of the entire aperture, but a workable 
compromise should be afforded with a hybrid system using the AMP technique to provide local control of 
a large number of small segments grouped on rigid "super segments” that would in turn be either 
interferometric or edge referenced to connect them and externally controlled only to compensate for the 
large-scale low-frequency flexing of the telescope structure and perhaps global tilt from turbulence In 
view of the limited number and low actuation frequency of the large segments, it might be useful to 
reassess the closed-loop multi-dither control approach from the early 70's, using target return and 
Iteration to maximize intensity. 


Program Recommendations 


Work on power beaming, at least in the near future, will have to be carefully prioritized in view of 
the amount that needs to be done and the limited funds available to support it. 

The FEL is one of the major elements of the system and considerable development will be 
needed to achieve short-wavelength operation at the conditions of interest, but that work, and especially 
the hardware needed to accomplish it, could prove to be expensive. It appears now that the pulse format 
of the induction FEL is incompatible with efficient photovoltaic power conversion and that the CW-like 
operation of the RF FEL will yield superior performance, but the latter may not as easily admit modular 
growth, and a new development plan with key milestones will have to be formulated. A program is in 
place for a lOOkW 10p RF FEL at Boeing and Los Alamos, Rocketdyne has most of the components for 
a IkW Ip system, research efforts are underway at Duke (with Russian support), Vanderbilt, and other 
laboratories, and there remains a wealth of induction FEL experience at LLNL, SRL, MIT and elsewhere 
A working group with knowledgeable participants from that community and from the government should 
assess the status of the RF FEL and develop realistic costs and schedules for a program leading to a 10 
MW 0.84p system. In addition, the potential cost and reliability advantages which made the induction 
FEL attractive in the past of course remain, and efforts should continue to be made to adapt that system 
to power beaming. 

The area that must receive the greatest attention in the near future is the adaptive optics, and 
the PAMELA hardware being acquired by MSFC for use there should provide an excellent test bed for 
sensors, algorithms, and actuators An important test of the concepts should be to mount that or similar 
equipment on a full-size truss telescope structure to demonstrate that compensation can indeed be 
accomplished on a system level, and the key technical results of that work should be the specification of 
component requirements and the levels of control ("woofer/midrange/tweeter") required for the most cost- 
effective operation. A working group with participants from at least MSFC, JPL, FORTH, KAMAN, AMP, 
and MITLL should periodically review that effort. An important area that has been largely neglected to 
date in this program is the synthetic beacon (or beacons) that will be required to provide the reference 
signal for the adaptive optics, and MITLL should be contracted to establish those requirements. 

The potentially high performance of the SRL/TRW electrodeless thruster could be extremely 
important for space transportation and work on it needs to be funded, but another source should be 
found to do so in view of its applicability to a broad range of NASA interests beyond power beaming 

Finally, high-level application studies should be conducted specifically to attract advocates from 
both the government and commercial communities, and an effort should especially be made to obtain the 
backing and resources of major aerospace and telecommunications contractors to whom power beaming 
could offer considerable long range benefits. 
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Attachments 


The attached technical reports were completed during or after the contractual period: 


A. 

B 

C 

D 

E. 

F 

G 


Optimum Segment Fitting to a Parabolic Surface 
Optimum Segment Fitting to a Parabolic Surface - Part II 
Optimum Gaussian Truncation 
EOTV Orbital Transfer 

Actuator Requirements for Turbulence Compensation 

Deterministic Control of Hexagonal Segmented Mirrors 
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OPTIMUM SEGMENT FITTING TO A PARABOLIC SURFACE 


Dr. Olenn W. Zeiders, The Sirius Group 


Calculations have shown that small segments with a fixed spherical figure can be used to 
approximate the surface of a parabolic mirror with satisfactory accuracy if the F/number is sufficiently 
large, and that there is very little benefit to using multiple spherical curvatures because most of the error 
is due to local astigmatism (different radial and azimuthal curvatures) of the parabola. However, 
advances in deterministic mirror figuring, such as LLNL's precision diamond machining, Itek's computer- 
controlled optical surfacing (CCOS), and Kodak's ion sputtering, may permit cost-effective, large-volume 
manufacturing of aspheric segments, and this paper explores the use of a limited number of different 
aspheric segments, carefully arranged to make best use of the basic hexagonal segment shape. The 
results indicate a significant benefit if the aspheric surfaces can be economically mass-produced . 


ANALYSIS 

A paraboloid is described by 

z = pr 2 =r 2 /4F, dz/dr = 2 p r = tan <p, 

and the segment surface coordinates ft, tj, Q are related to fixed-space coordinates (x, y, z) for a 
segment defined to be located on the x-axis by 

x = r + % cos <p - £ sin <p , y = t| , z = p r 2 + tjsin<p+£ cos q>, 

from which it follows that £ = p ft 2 cos 2 q> + q 2 ) cos q> = P ft 2 + if 2 / Y 3 . 

neglecting terms of order d/F relative to unity and where 

y 2 = sec 2 <p = 1 + (2 p r) 2 . 



It can be seen from inspection of these equations that the local radial and azimuthal radii of curvature 
are given by 2 p R r = y 3 and 2 p R e = y. 

Now, letting the segment surface figure be astigmatic with radii R x and Ry defined in cartesian 
space, the mean square distortion across a circular segment of diameter d (d = 1 .(KOI x the face-to-face 
width of a hexagon of the same area) is given by 
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s 2 = (4/nd2)J^ ,2 | [5 0 + p2 { cos 2 (0 + <p) / 2R X + sin 2 (0 + <p) / 2Ry - p cos 2 ^ / y 3 - p sirAp / y }]2 p d© dp 


= 5 0 2 + 8 0 d 2 {1 / 2R X + 1 / 2R y - (y 2 + 1) / y 3 ) / 8 + d 4 [(1 / 2R X - p / y 3 ) 2 + (1 / 2R X - p / y3)(1 / 2R y - 
1/2R x )-(1/2R x -p/y 3 )P( Y 2 -1)/y 3 +3(1 /2R y -1/2R x ) 2 /8 + 3p 2 (y 2 -1) 2 /&y 6 - 
(1 / 2 R y - 1 / 2Rx) p (3 - 2 sin 2 ©) (y 2 - 1) / 4y3. 

8 0 is a simple “piston" correction that can be locally applied, and the optimum value is 

8 0 = -d 2 {1/2R x + 1/2R y -P(y 2 +1)/y 3 )/16, 

sothat 8 2 = (d 4 /192)[(1/2R x -p/y 3 )(1/2Ry-p/y) + 3{1/2R y -1/2R x -p(y 2 -1)/y 3 } 2 /4 

+ 2 p (1 /2R y -1 /2R x )sin 2 0(y 2 -1)/y 3 . 

For the special case of spherical segments, the two radii are equal, and the equation reduces to 

8 2 = d 4 [{1 / R - p (y 2 + 1) / y 3 ) 2 + 2 p 2 (y 2 - I) 2 / y®j / 768. 

Note that optimal variation of the spherical radius of curvature, i.e. pR = y 3 /(y 2 +1), would still result in 
an appreciable astigmatic contribution from the second term. Averaging over a region from D-j/2 = a D/2 
to D/2, the overall mean square variation of the output wavefront is therefore given by 

d/2 0 

< Sw^ > = [4 / {* D 2 (1 - a 2 )}] | 8 2 4 serAp 2* r dr = (8 / p 2 D 2 ) J 8 2 y 3 cfy 

0 0i 

= (d 4 / 96 D 2 ) [(y 2 - 1 ) / (4 p 2 R 2 ) - 2 {(y 3 - 1 ) / 3 + y - 1} / (P R) + (y 2 - 1 ) / 2 + 2 ln(y) + 

(Y 2 - 1) / (2 y 2 ) + 2 {(y 2 - 1) / 2 - 2 ln(y) + (y 2 - 1) / (2 y 2 )}], 

from which it follows that the optimum constant spherical radius of curvature is given by 

(y 4 - yi 4 ) / (4 P R) = (^-ri^/a + y-y!, 

sothat < Syy 2 > = (P 2 d 4 /96)[3{1 + 1/(y 2 y 1 2 )}/2-2ln(y)/(y 2 -y 1 2 )- 

4 (y 2 + m + Yi 2 + 3) 2 / {9 (y 2 + Yl^y + Yi) 2 )] 

= (p 2 d 4 /192)[(1 +a 4 )e 4 - (19 + 4a 2 + 7a 4 + 18a 6 )e 6 /12+ ... ] 
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where e = 0 D = D / 4F. The segment size requirements for a full aperture (a = 0) are shown in Fig. 1 for 
both a continuously variable optimum spherical figure and a constant optimum spherical curvature. The 
two differ by only about 10% because of the dominance of the astigmatic term, so it would not be cost- 
effective to vary the spherical curvature either continuously or piecewise across the aperture. 

Returning to the case of an astigmatic correction and integrating that expression over the 
annular region, we obtain 


<5w 2> = [d 4 / {48 D 2 (1 - a 2 )}] [(y 4 - yf 4 ) / (8 P 2 R x Ry) - (y 3 - Tl 3 )/ <3 p R„) - (y- Yl ) / (P R y ) + 
3 (/-Yi 4 ) {1 / (2 p R y ) - 1 / (2 p Rx)} 2 /8 + 3 (y 2 -^ 2 ) {1 + 1 / (y 2 n 2 )} /4 - ln(y/ yl) - 
{(r 3 - n 3 ) / 3 - y + yi > {1 / (2 p Ry) - 1 / (2 p RJ} {1 + 2 (sin 0) / 6}], 


from which the optimum curvatures are found to be given by 


(Y 4 - ri 4 ) I (4 P Ry) = (y 3 - y- 1 3 ) / 3 + y - y<| + {(y 3 - y-j 3 ) / 3 - y + y-j} (sin 9) / 6 


(y 4 - Yl 4 ) / (4 P Rx) = (Y 3 - Yl 3 ) / 3 + y - y-j - {(y 3 - y 1 3 ) / 3 - y + y-|} (sin 9) / 9. 

<5 w 2 > = (p 2 d 4 /96)[3{1 + 1/(y 2 y 1 2 )}/2-2ln(y)/(y 2 -y 1 2 )- 
4 (y 2 + YY1 + Yl 2 + 3) 2 / {9 (y 2 + yi 2 )(y + y^) 2 }} - 
8 (Y 2 + yyi + y-| 2 - 3) 2 {(sin 9) / e} 2 / {9 (y 2 + y 1 2 )(y + y-,)^] 

= (p 2 d 4 / 192) [(1 + a 4 ) e 4 - (19 + 4 a 2 + 7 a 4 + 18 a 6 ) s 6 / 12 + ... 

- (1 + a 2 ) {(1 + a 2 ) e 4 / 2 -((1 - a 2 ) 2 / 6 + (1 + a 2 ) 2 ] e 6 / 2 + ...}{(sin 0) / G} 2 ] 


This is plotted in Fig. 2, again as segment size vs. f/number for 1/20 wave distortion Q 0.8 p. with a 12 m 
aperture, and the results are shown for astigmatic segments with 30° and 60° included angles and for 
annuli with a = 0.5, 0.8, and 1 ; the equivalent results for segments with constant spherical figure over the 
annuli are shown for comparison. Unlike the all-spherical case where little is gained by segmentation, 
relatively coarse segmentation with astigmatic elements can offer significant benefit, but it should be 
noted that very fine annuli (i.e., a -* 1) offer a benefit only with very fine angular segmentation as well (9 
0), and that would not be likely to be cost-effective. 

Key to mass-production of aspheric elements by any of the new deterministic surface figuring 
techniques is the amount of material to be removed after an initial spherical figure is applied. It follows 
that the maximum aspheric thickness to be removed with the optimum astigmatic segments is 

d 2 (1 / Ry - 1 / Rx) / 8 = (d 2 / 12F) (y 2 + yy-| +y-, 2 - 3) {(sin 9) / 6} / {(y 2 + y-^y + y-j)} 

= (e 2 d 2 / 16F) (1 + a 2 ) {sin 9/0} + ... 


and 

Thus, 
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As an example, 0.54 p would be the maximum aspheric thickness that would have to be removed from 4 
cm segments for a 12 m, f/1 mirror having an astigmatic zone with 0 = 00° and a = 0.5. While large 
compared to the 0.02 p figuring tolerance at 0.8 p wavelength, it's small compared to the d 2 /16F = 8.3 p 
spherical well depth, and it should be possible to accurately remove the material in a single pass. 

SYSTEM OPTIMIZATION AND RESULTS 


A set of computer experiments were conducted to determine the "optimum" configurations for 
use with hexagonal mirror segments, noting that 60° sections are repeated, and that integer fractions of 
60° simply require that the primary axes be rotated relative to an edge. Every possible combination for 
up to six different individual segment configurations were tried, and the results for the best solutions are 
given in Table I for an annular primarily mirror with an inner diameter of 10% of the outer and, except for 
the spherical reference case, for an inner astigmatic ring with identical segments (ie., repeated every 
80°.) The results are again presented in terms of the required equivalent circular diameter for 0.04 p 
RMS wavefront distortion over a 12 meter mirror (in bold type), and the inner diameters of the astigmatic 
zones are given in italic type. The geometries are shown in Fig. 3. 


Number of 

Number of 

Angular 



F/Number 



Different 

Different 

Size of 






Figures 

Segments 

Outer 

0.5 

1.0 

1.5 

2.0 

2.5 



Zones 






1 

i 

Spherical 

0.7979 

2.1136 

3.8317 

5.8715 

8.1875 

1 

i 

60 

0.8995 

2.3608 

4.2720 

6.5417 

9.1193 

2 

2 

60 

0.686 

0.704 

0.708 

0.700 

0.710 




1.0625 

2.8077 

5.0873 

7.7939 

10.6673 

2 

3 

30 

0.663 

0.683 

0.688 

0.680 

0.600 




1.2204 

3.2036 

5.7968 

8.8765 

12.3740 

3 

4 

60 

0.512 

0.534 

0.530 

0.540 

0.540 



30 

0.732 

0.782 

0.757 

0.758 

0.750 




1.2997 

3.4237 

6.1996 

9.4956 

13.2388 

3 

5 

30 

0.532 

0.554 

0.558 

0.560 

0.560 



30 

0.700 

0.805 

0.800 

0.810 

0.810 




1.3769 

3.6271 

6.5676 

10.0594 

14.0246 

3 

6 

30 

0.522 

0.543 

0.540 

0.550 

0.550 



20 

0.773 

0.700 

0.704 

0.705 

0.705 




1.4630 

3.8406 

6.9495 

10.6416 

14.8345 


Table I Optimum Configurations with Astigmatic Segments 

At a given f/# - and if other considerations should allow it - the segment size could be increased by 
81% by using 6 different astigmatic segments (but with only 3 different figures), thereby decreasing the 
total number of segments by a factor of 3.3. Alternatively, were a system originally to be configured at 
f/1. 5 on the basis of spherical segments, a similar design change would permit the use of a f/1 .0 
system, thereby decreasing the height of the secondary mirror structure by 33%. Furthermore, the inner 
zone probably doesn’t require aspherics, because it has been found that the use of spherical segments 
there would only increase the minimum allowable segment size by a few per cent at most 
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Note that results have been presented for a particular set of design conditions (0.04 p RMS 
wavefront distortion over a 12 meter mirror), but the maximum allowable equivalent diameter can be 
scaled to other conditions by the fact that it is proportional to (SyyD)^ for all other design parameters 
remaining unchanged. The dependence upon allowable distortion can be especially important because 
the useful power on target scales as exp[-(2*&fl,)2] when the wavefront error correlation lengths for the 
disturbances are small (as they are here). Thus, the factor in the exponential loss scales as d* for this 
geometrical effect, whereas it scales only as for Kolmogorov turbulence, so it produces a much 
more dramatic loss than does atmospheric turbulence as the segment diameter increases beyond a 
certain value - and that value may be increased significantly by the use of aspherics. 


Although quite attractive from the performance standpoint, the ultimate cost effectiveness of 
aspherics wilt depend upon the ability to produce them economically in large quantities. 


The formula manipulations in this paper are straightforward but painfully lengthy, and details have been 
presented only to indicate the approach or to give important results. Details are available from the 
author upon request. 
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SEGMENT DIAMETER d, CM. SEGMENT DIAMETER d, CM. 
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Fig. 1 Spherical-Figured Segment Requirements 
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Fig. 2 Astigmatic Segment Requirements 
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TVvo Aspheric Figures, 
Two Orientations in 
Outer Ring 


Three Aspheric Figures, Three Aspheric Figures, 
Two Orientations in Two Orientations in Eacn 

Outer Ring of Two Outer Rings 



Three Aspheric Figures, 
Three Orientations in Outer 
Ring and Two in Middle 


Fig. 3 Calculated Configurations 
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OPTIMUM SEGMENT FITTING TO A PARABOLIC SURFACE - PART II 

Dr. Glenn W. Zetders, The Sirius Group 


A subject that continues to reappear in the SELENE program is the question of fitting hexagonal 
segments to a parabolic surface. This was briefly addressed at the 2/92 JPL Wavefront Sensing and 
Control Workshop by the author and in an analysis by D. Korsch who showed that non-interference 
could be assured by simply projecting and tilting the segments from the reference plane onto the 
paraboloid. However, because radial arcs lengthen while azimuthal ones remain unchanged under such 
a transformation, that simple fitting procedure produces increasingly large radial gaps away from the 
center while leaving essentially no azimuthal ones, and that can lead to unnecessarily large energy loss 
and, perhaps, difficulties in edge sensing. Here we begin with the same proposition, i.e. projecting the 
segments from the reference plane onto the curved surface, but we permit outward radial shifts to relieve 
the azimuthal "crunch", we specifically require that the gaps be locally symmetric (i.e., azimuthal = 
radial), and we consider the benefits of radial zones with slightly different-sized segments to close the 
gaps. 


If we allow segments to be shifted outward by a distance d(t) and to have an effective diameter 
d(r), the mean azimuthal gap for segments initially at a distance rfrom the center is 


Ae=[2n(r + 5)-Nd]/N, Nd 0 = 2?rr => Aq = (1 + 8 /r) d 0 - d 
and the radial gap between segments initially at r and those at r + d Q is 

A p = (r + 8 - d« cos 6)^ do - (r + 5 + tV2 cos 0) r = d 0 + do - d 0 5 / r - d 0 /2 (d cos 0) 


where tan 0 = r / 2F, and where we have modeled the discrete segments by a continuous distribution, 
valid in the limit of very small segment size. Then, we find 


Ap - Aq = d 0 °®/0 r - d Q 8 / r + d/2 (r / 2F) 2 - do 12 + (small terms ) , 

and we will henceforth let A = Ap = Ag to maintain local symmetry. 

Case I: A - 0 Letting 8 = 0 @ r = D/2, d/do = S/r+ 1 = expKE^M-r^/lOF 2 ] 

which would require a continuous distribution of segment sizes for which the central segment would be 
1 .6% larger than those at the periphery for f/1 .0. 


Case II: d ■ d 0 Again with 5 = 0 @ r = D/2, 8 / r = A / d = (d 2/4 - r 2 ) / 16F 2 . 

This gives a maximum gap Zl = d (D / 8F) 2 Q r = 0, which is 625ji (0.025") for d = 4 cm and f/1 .0, but 

this would reduce to 31 2p if 4.031 cm segments were used within r = D/V2, and the gap would decrease 
further with increased zonal segmentation. 
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OPTIMUM GAUSSIAN TRUNCATION 


O. W. Zeidars, Tha Sirius Group 


The output beam from a free -electron laser tends to be a freely-propagating one with virtually no 
edge effects of edge-clipping, and that leads to an essentially infinite Gaussian profile which must be 
processed by the finite-sized transfer optics and the BTOS for final propagation to the target Excessive 
power will be wasted if these elements truncate the beam at too small a radius, while excessive 
diffraction will occur if the truncation radius is too large, so an optimum ratio of beam/waist size occurs 
for which the peak far-field intensity (i.e., ideal Strehl ratio) is maximized. 

The profile of an infinite gaussian beam is described by 


I = l 0 expff/a) 2 /?} sothat = 2*jlrdr = 2n a 2 l Q . 

o 


If the beam is truncated at a radius R, the peak far-field intensity is then proportional to 


R R 2 

[2 1 r Vl dr J 2 = [ V(P 00 /2sa 2 ) J exp{-(r/2a) 2 } dr 2 ) 2 

R8 (RS) 2 


where RS is the radius of the secondary mirror. Then, optimizing with respect to the beam/waist size 
ratio, we have (symbolically) 


Y = [(exptfSx) 2 } - expfx 2 }) / x ] 2 (x = R / 2a) 


sothat 


dy / dx = (2 x 2 ( -6 2 expf (Sx) 2 } + expfx 2 }) - (expffSx) 2 } - expfx 2 })] 2 Vy / x 2 , 


is zero for optimum conditions, and this requires that 


(1 + 2 x 2 ) expfx 2 } = (1 + 2 8 2 x 2 ) expKSx) 2 }. 


The resulting optimum beam size r/a, the Strehl ratio relative to a full (S = 0) infinite beam with 
the same waist size, and the fractional powers 1 - expf (5R/a) 2 /2} within the inner radius, expf (R/a ) 2 /2} 
beyond the outer [also the edge/center intensity ratio], and the net are shown in the table below as 
functions of 5: 


d 

R/a 

Strehl 

Inner Power 
Outer 
Power 
Net Power 


0 

0.05 

0.10 

0.15 

0.20 

0.25 

2.2418 

2.2355 

2.2166 

2.1870 

2.1490 

2.1046 

0.71533 

0.71020 

0.69501 

0.67098 

0.63967 

0.57463 

0 

0.00623 

0.02427 

0.05230 

0.08823 

0.12926 

0.08104 

0.08219 

0.08572 

0.09149 

0.09935 

0.10919 

0.91896 

0.91158 

0.89001 

0.85612 

0.81242 

0.76155 


This corresponds to truncation slightly beyond the e ~ 2 intensity points, and it might be expected that an 
optimization considering the cost of the telescope would result in a somewhat smaller value of R/a. 
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EOTV ORBITAL TRANSFER 


Dr. Glenn W. Zeiders, The Sirius Group 


In response to questions concerning the long Earth-orbital transfer times being reported, this report 
reviews some of the key principles that are important to the power beaming program. In the first part, the 
trajectory parameters for two-phase constant-acceleration thrusting from low orbit to a higher one are calculated 
first tangential-directed for rapid energy addition, then normal-directed (constant energy) to meet the angular 
momentum requirement The model is not meant to represent optimum conditions, but it does simulate an 
interesting flight profile, clearly shows the most important behavior, and gives some useful systems-level results: 

1 Continuous thrusting is far more costly than a Hohmann two-impulse, half-orbit transfer 
from the AV standpoint. 

2. Acceleration is probably the most important single factor determining the mission 
profile. The required AV increases only slightly with acceleration level, but mission duration 
decreases dramatically High thrust and/or low weight are suggested to safely transit the Van 
Allen belts. 

3. At very low acceleration levels, the vehicle spirals slowly outward about the central 
body, and the radial velocity is sufficiently low that very little final maneuvering is required for 
injection into the new orbit. The overall AV can be shown to approach the difference of the two 
circular velocities, and this leads to a particularly simple closed-form solution which should prove 
useful for overall systems analysis as well as for interpretation of more detailed computer 
simulations. 

In view of the preceding, the second part examines the requirements on the system specific weight and the 
specific impulse, and concludes that it is important that the former be reduced as much as possible whereas the 
latter should be chosen judiciously with flight time in mind. High specific impulse saves propellant mass at the 
expense of thrust (acceleration), and, when the propellant mass becomes a small fraction of the total package (as 
it does at high l s ), the specific impulse may be able to be reduced significantly to dramatically reduce flight time 
without appreciably impacting overall system mass. Calculations for realistic conditions yield trip times 
considerably shorter than reported by Boeing. 


1. Continuous Thrusting Transfer 

The planar equations of motion in polar coordinates and 
dimensionless form are 

d 2 r / dt 2 - v^/r+1/r^ = a cos y 
d(rn)/dt = rasiny 

where r = r/r 0 t = t^/r 0 3 

v = V sin q> y r 0 / n a = A r 0 2 / ji 



/ 

l _ 

\ / 


Fig. 1 Maneuver Geometry 
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and where the Earth gravitational constant p = 3.9853 10 5 km 3 / sec 2 . For transfer from low-Earth orbit (LEO) 
at 500 km attitude, the initial radius r 0 is 6878 km, and the other normalizing factors are -f 0 3 /p = 904 sec = 0.251 
hrs for time, ^p/r 0 = 7.61 km/sec for velocity, and p/r 0 2 = 8.42 m/sec 2 (0.860 G's) for acceleration. 




®T 

Tt 

t T 

dT 

a N 

tfsl 

dN 

tT + tN 

dT + 

®T^T + ®N*N 

— 










dN 



6 

.001 

5.95 

591 

243 

.0081 

3.59 

.258 

595 

243 

.620 



.003 

5.69 

196 

81.2 

.0088 

8.57 

.604 

205 

81.8 

.663 



.01 

4.83 

57.2 

24.6 

.0126 

13.8 

1.06 

71.0 

25.6 

.746 



.03 

3.45 

17.6 

8.55 

.0194 

16.3 

1.48 

33.9 

10.0 

.845 



.1 

2.11 

4.18 

3.16 

.0314 

15.5 

1.67 

19.7 

4.83 

.904 



.3 

1.18 

1.21 

1.31 

.0347 

16.3 

249 

17.6 

3.80 

.931 



1 

1.02 

.355 

414 

.0352 

16.7 

3.03 

17.1 

3.45 

.942 









Hohmann: 

20.6 

3.14 

.499 


5 

.001 

4.97 

553 

240 

.0087 

2.56 

.238 

555 

240 

.575 

— 


.003 

4.85 

183 

80.1 

.0128 

4.81 

.441 

188 

80.5 

.612 



.01 

4.32 

53.9 

24.2 

.0163 

9.30 

.904 

63.2 

25.1 

.691 



.03 

3.21 

16.7 

8.40 

.0242 

11.3 

1.40 

29.1 

9.79 

.780 



.10 

2.01 

4.00 

3.10 

.0410 

11.8 

1.60 

15.8 

4.70 

.883 



.30 

1.16 

1.17 

1.27 

.0450 

13.0 

2.49 

14.1 

3.76 

.933 



1.0 

1.02 

.342 

.398 

.0455 

13.3 

3.01 

13.7 

3.41 

.950 









Hohmann: 

16.32 

3.14 

.480 


4 

.001 

3.99 

500 

234 

.0110 

1.49 

.201 

501 

235 

.516 



.003 

3.93 

166 

78.2 

.0168 

2.82 

.360 

169 

78.6 

.546 



.01 

3.68 

49.1 

23.6 

.0244 

5.18 

.685 

54.3 

24,3 

.617 

— 


.03 

2.68 

15.5 

8.15 

.0307 

8.67 

1.30 

24.2 

9.44 

.730 



.10 

1.88 

3.72 

2.99 

.0559 

8.60 

1.55 

12.3 

4.53 

.852 



.30 

1.14 

1.11 

1.21 

.0600 

9.89 

2.47 

11.0 

3.68 

.926 



1.0 

1.01 

.326 

.377 

.0607 

10.3 

2.99 

10.6 

3.37 

.950 









Hohmann: 

12.42 

3.14 

.449 


3 

.001 

3.00 

423 

222 

.0131 

.849 

.018 

423 

222 

.434 



.003 

2.99 

141 

74.1 

.0745 

.410 

.080 

141 

74.2 

.453 



.01 

2.88 

41.8 

22.4 

.0406 

2.45 

.496 

44.2 

22.9 

.518 



.03 

245 

13.4 

765 

.0391 

5.45 

1.18 

18.9 

8.83 

.616 



.10 

1.69 

3.30 

2.80 

.0800 

5.73 

1.47 

9.03 

4.27 

.788 



.30 

1.10 

.995 

1.09 

.0829 

7.08 

2.48 

8.08 

3.57 

.885 



1.0 

1.01 

.291 

.332 

.0839 

7.46 

2.98 

7.75 

3.31 

.918 









Hohmann: 

8.89 

3.14 

.393 

— 

2 

.01 

1.98 

29.2 

18.8 

.0445 

1.01 

.370 

30.2 

19.2 

.337 



.03 

1.90 

9.69 

6.41 

.0403 

2.62 

.983 

12.3 

7.39 

.396 



.10 

1.35 

2.42 

2.29 

.1153 

3.45 

1.53 

5.87 

3.82 

.640 



.30 

1.05 

.762 

.830 

.1103 

4.66 

2.56 

542 

3.39 

.772 



1.0 

1.00 

.225 

.250 

.1113 

5.00 

2.98 

5.23 

3.23 

.782 









Hohmann: 

5.77 

3.14 

.284 


TABLE I CONTINUOUS THRUSTING ORBITAL TRANSFER 


In-plane maneuvers are accomplished by appropriately programming the magnitude A and direction y of 
the acceleration, and this analysis specifically assumes constant tangential acceleration from a low circular orbit 
to the required final energy for a higher circular orbit, followed by constant normal acceleration to the required 
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angular momentum for the final orbit. In the first operation, the energy increases linearly with path length 
according to 


E = V 2 12 - 1 / P = <x T S - 1/2 , S = I V dt , AS = (1 - r 0 / r F ) /2 a T , 


and the trajectory is described by 


p d2p / dS 2 = (1 - 1 / p V 2 ) [1 - (dp / dS) 2 ] where p = 1, V = 1, and dp / dS = 0 @ S= 0. 
The trajectory for the second maneuver is described by 


p d 2 p / <h 2 = [1/ p - 1/ pp - (dp / dr) 2 ] - a^| ^ [1 - (dp / dr ) 2 1 (21 p - 1/ pp)] 


where the acceleration required to match the boundary conditions can be shown to be given by 


VPp - Hj = 3a N pp 5/2 (ji/ 4 - 2/3 - sin'^pj^pp + (1 + 2/3 pj/2pp } >/pj/2pp ^[1 - pj/2pp]) 


where Hj = p V ^[1 - (dp/dS) 2 ] is the angular momentum at the end of the tangential-thrusting phase. 

These equations have been numerically integrated on a PC, and the results are presented in Table I 
where the last three columns represent, respectively, the total normalized transit time, the total angular traverse 
(radians) about the central body, and the total normalized impulse. These quantities have in turn been plotted in 
Figure 2, and the corresponding dimensioned values are shown for a 500 km lower orbit (note that pp/p 0 = 6.1 for 
a flight from LEO to GEO.) Also shown are the equivalent values for an ideal Hohmann transfer; they become 
substantially more favorable as the altitude change increases, but the maneuver requires infinite-acceleration 


Total Impulse, km/sec 



Flg.2 Orbital Transfer Requirements 
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impulses at the beginning and end of a half-orbit and is not suitable for electric propulsion alone. The oscillations 
that may be noted in the results are probably not calculational artifacts, but rather are most likely due to the 
tendency of the vehicle to "hunt' about the mean trajectory; such behavior is clearly evident from the presence of 
the trigonometric functions in Benne/s power series solution 


p - 1 + 2 a (S - sin S) + 4 a 2 [S 2 + S/2 (S cos S - sin S) + 2 (cos S - 1)] + . 


for the case of constant tangential acceleration from a circular orbit. 

It should be obvious from Figure 2 that the very long-duration missions being projected by Boeing reflect 
acceleration levels well below 0.001 and correspond to operation above the upper edge of the graph. For those 
conditions, the trajectory approaches a gradual spiral with many orbits around the central body, and simple 
closed-form results can be obtained in that case by noting that the velocity is essentially tangential and equal to 
the local circular value, l.e. pV 2 s 1 . Thus, 


2aSs 1-1/p and x = jdS/V s ( 1 - 1/Vp)/a 


so that the velocity increment is just the difference between the initial and final orbital velocities the vehicle spirals 
outward with 

p s 1 / ( 1 - a x ) 2 , 


and the mission duration for a given orbit change is essentially inversely proportional to the acceleration. 


These results should prove quite useful if one wishes to 
quantitatively examine the effects of laser propagation loss and/or 
intermittent thrusting as the vehicle passes through "irradiation cones." The 
geometry for the latter is shown in Figure 3, from which it follows that the 
exposure fraction at a particular altitude is given by 


<P / * = [6 - sin'^sin e / (1 + h / R 0 ) } ] / * 


for each ground station; this is plotted in Figure 4. The mean acceleration is 
proportional to the value, so the time spent near that altitude is inversely 
proportional to it. For a maximum angle of 60° from zenith, a vehicle at 500 Fig. 3 Irradiation Geometry 
km altitude would be accessible from each ground station 3.65% of the time 

(6.33% at 75°!), and it would rise to 6.40% at 1000 km and 10.42% at 2000 km. Although this matter will tend to 
heal itself as the program progresses and more ground stations become available, the real problem with low- 
acceleration transfer from LEO will remain: low attitudes where the radial velocity is the lowest (and the "masking" 
is the greatest) can also be the most lethal from the standpoint of the Van Allen radiation belts in which charged 
particles are trapped by the Earth's magnetic field. The inner belt, extending from the surface to about 2R Q and 
peaking at about 1.5R 0 (3000-4000 km altitude), is composed primarily of penetrating protons of cosmic origin, 
whereas the outer belt, which extends nearly to GEO, is made up of low-energy electrons. A vehicle rising from 
LEO with constant tangential acceleration would spend at least 30% of its flight time in the more-dangerous lower 
radiation belt, and this certainly suggests the need for a high-thrust system to minimize the time spent at low 
altitudes, either using an initial kick from an auxiliary rocket or using a laser-propelled vehicle with power from 
relay satellites rather than from the surface and, perhaps, with degraded specific impulse to Increase the thrust. 



Acceleration is the key parameter determining mission duration, and 
short flight times will therefore require high thrust and/or low mass. 
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Ag.4 Coverage per Ground Station 


2. Vehicle Considerations 


The parameters of the vehicle are related by 

F = m l s g 0 and P B = mv 2 /2 = m(l s g 0 ) 2 /2 = Fl s g 0 /2 

where F = thrust l s = specific impulse 

m = mass flow rate P B = power in the thrust beam = q P. 

Thus, dv / dt = F / M = m dv / dm = F / (l s g Q ) dv / dm, 

which Integrates to the familiar rocket equation 

AV = l s g Q ln( M 0 / M) = l s g 0 /n[1 / (1 - M prop / M 0 ) ] 

for constant specific impulse. Then, assigning a fixed mass Mpix (payload, etc.) and individual time-dependent 
(shielding, etc.) and power-dependent masses and M POWER' the inltial mass is 9 iven b y 

[M F |x + (MpOWER /p ) PI > M o = 1 * {1 * exp[-(AV/l s g 0 )IK1 + {(Mtime/^) (AV) 2 /(2Pti)} (l s go/AV) 2 ] 


and the mission duration t by 

[M F |x + (MPOWER^) PI (AV) 2 /2 Pt, / x = (AV/l s g 0 ) 2 / (exp[(AV/l sgo )] - 1} - {(M TIME /x) (AV) 2 /(2Pq)}. 
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These results are plotted in Figure 5, from which it can be seen that the transfer time is minimized for Lg^/AV = 
0.6275 (I.e., I s < 300 seconds for all LEO/GEO transfers.) Although the curves for various values of the time- 
dependent mass parameter are interesting, their utility is limited because the parameter tends to be small; as a 
point of reference, using Grant Logan's shielding calculations presented at the 2/91 Beam Power Workshop at 
LeRC (2.5 tons for a 10 day transfer with 3.7 tons of propellant, l s = 2500 sec, and AV = 4.3 km/sec), we find the 
parameter to be only 0.02. Nevertheless, the analysis Is important because it shows that the gross weight (often 
used as a measure of overall performance for a space vehicle), varies weakly for large values of the specific 
impulse whereas the mission time varies essentially linearly with l s in that region, so virtually any weighting of 
mission duration is likely to drive one to I^q/AV =2-3 for best operating conditions. To put Figure 5 in 
perspective with real parameters, consider the following: 


Then, 


l s = 1200 sec. 


AV = 4.536 km/sec (500 km to GEO at low a) 

P = 500 kW with i] = 0.5 

•s9</4V = 2.593 M powER = 2 5 tons 


MpowER^ 3 = ® kgm/kW 
Mp|x = 2.5 tons 

WTIME = 0 

Mrrop = 2.35 tons 


M 0 = 7.35 tons 


At = 7.53 days 


F = 42.5 newtonsF/M 0 g 0 = 0.0006 


This is about twice as long as the minimum time - which would require substantially more propellant - and is only 
47% heavier than would be achieved with Infinite l s - and infinite mission time. 



FlQ- 6 Effects of Specific Impulse on Mission Duration and Goss Weight 

Although not specifically plotted, one should note the effect of specific power M' = M P owER /p in tt 16 
equations. It always contributes negatiely by increasing the overall system mass, and it serves to limit the useful 
power to Mpix/M’. For a fixed mass Mpix of 2.5 tons, a value of M' = 5 kgm/kW would effectively limit the useful 
power to 500 kW, and 1 kgm/kW would be needed for 2.5 MW. 


Increased specific impulse reduces propellant mass at the expense of thrust 
(acceleration.) Once the propellant mass becomes a small fraction of the total 
system, specific impulse may be reduced substantially to decrease mission time. 

A lower value of specific power may be required to efficiently use the high 
input power available with laser beaming. 
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ACTUATOR REQUIREMENTS FOR TURBULENCE COMPENSATION 


Dr. Glenn W. Zeiders, The Sirius Group 


A simple methodology, based on the work of Noll (Ref. 1) and Greenwood and Fried (Ref. 2), is presented 
for estimating the actuator requirements for adaptive optics correction of global tilt and residual piston error arising 
from Kolmogoroff atmospheric turbulence, and predictions are given for SELENE. The results of G&F for the 
frequency characteristics are specifically applied to ground-to-space conditions for which the segment sizes are very 
small compared to the full aperture, and results are given for several Important crosswind distributions; It is found that 
the velocity profile has a strong effect on the power spectra and that high slew rates in particular significantly 
Increase the required high-frequency response and accentuate the effects of high-altitude turbulence. 


Atmospheric turbulence disturbs the air density and index of refraction along an optical path, large turbules 
primarily causing low-frequency centroid wander of the central diffractive lobe while small ones reduce the peak 
intensity, scattering energy over wide angles. The amount of wander and scatter for a uniform plane wave traversing 
a path along z through full-range Kolmogoroff turbulence are characterized by the scale 


r 0 = {0.423 k 2 \ C n 2 (z) dz}‘ 3/5 


where k = 2nfk is the wave number of the radiation; for the Hufnagel-Valley turbulence profile (Fig. 2) and a path 
inclined at an angle 8 from zenith, rg = 11.1 {A^ 2 cos 8} 3 ^ cm where 73% of the contribution arises from the "ground 

layer" and only 9% is due to the "tropospheric spike.” It is now well-known that appreciable improvement can be 
realized in both imaging and beaming applications if the disturbed wavefront is properly corrected over dimensions of 
r 0 or less. As long as the remaining distortion is not too large, the relative peak intensity (Strehl ratio) is reduced to 
exp(-<o 2 >) where <o 2 > is the variance of the corrected wavefront from its mean. 

Again for full-range Kolmogoroff turbulence with a plane wave and further assuming a uniformly-illuminated 
circular aperture of diameter D, Noll (Ref. 1) exactly evaluated the residual variances after progressive Zemike 
polynomials were removed, and showed that the mean variance for the "uncorrected" wavefront is 


Piston-Corrected: <o 2 >p = 1 .030 (D/rg)®^ 


after time-dependent piston is taken out to give zero mean (unnecessary In practice, but convenient for calculations), 
while the residual variance of a fully tilt-corrected wavefront is 


Piston/Tilt-Corrected: <o 2 >p/j = 0.134 (D/fg)^. 


The tilt is equally probable in all directions, so the total time-averaged variance can be represented symbolically by 6 2 
= a 2 + eV, and the spatially-averaged variance is then <S 2 > = (4/jrD 2 ) / 2«r JS 2 dr = a 2 + ^0^/8. Thus, the 
variance of the tilt alone is 


<0 2 > = 1.792 (D/r 0 ) 5/3 QJnD) 2 
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The aperture diameter in the above formulas is actually quite arbitrary, and the results can be applied 
equally well to the full A/O problem by using the full aperture size D to calculate <&■> for the global tilt corrector and 
the resulting wavefront variance <cP->p/j which is to be corrected by piston motion of the subapertures, while the tilt 
correction and resulting variance for a subaperture can be calculated with the same formulas by using its diameter d 
instead. Thus, 


Tilt Mirror Requirement. 
Subaperture Requirement: 
Residual Wavefront Error: 


<e2> = 0 182 (D/r 0 ) 5/3 (X/D) 2 
<o 2 > p/t = 0.134 (D/r 0 ) 5/3 
<o 2 > P/T = 0 134 (d/r 0 ) 5/3 


<02> = 0.182 (d/r 0 ) 5/3 (X/d) 2 


so, in terms of allowable residual RMS wavefront error, the required subaperture size would be 


Srms : 

X/5 

X/ 10 

XI 20 

X/40 

Max. Strehl: 

0.21 

0.67 

0.91 

0.98 

(^P/T-CORRECTED 

4.39 

1.91 

0.83 

0.36. 


The primary advantage of global tilt removal is that large control displacements are required only at 
relatively low frequencies, whereas, were it not done, the high-frequency controls for the individual 
segments would have to be capable of mean strokes of ( 1. 030/0. 134 = 2.77 rimes that for residual 
correction alone, and the edge segments would have to move further yet. That may or may not be a 
problem; the cost/benefits of the separate global tilt correction system would depend upon the particular 
implementation, and would require a separate assessment based specifically upon it. 

An ideal gaussian distribution is unbounded (although real effects will produce real limits in practice), so the 
translation from variance to full stroke tends to be a matter of individual choice, but the minimum realistic value to use 
as a base would appear to be V8, the reason being that it is the appropriate factor for a sine variation, i.e 


2it 

< ♦ > 2 = 1 ( 0/2 sin e ) 2 d0 / 2« - O 2 / 8, 

0 

and this should next be multiplied by "safety factor of perhaps 1.5 - 2 to avoid "going up against the stops", then 
divided by another factor of two to account for the fact that phase is shifted by twice the displacement of a mirror 
upon reflection; thus, the overall multiplier would be 2.1 - 2.8 to go from root variance to P-P stroke (Lincoln Lab has 
consistently used 2.5 in the past, and it would seem appropriate to adopt that for SELENE as well.) 

For a typical SELENE scenario for which the wavelength would be about 0.8p to optimize the 
response from GaAs photovoltaic cells, the coherence scale size tq with Hufnagel-Valley turbulence would 
be 8.5 cm at zenith, but that would decrease to only 3.8 cm at 75P. The larger number of subapertures 
required for correction at large zenith angles could prove to be a very costly penalty to pay for increased 
angular coverage, and it is likely that a proper optimization might well call for higher laser power instead to 
provide the same delivery at reduced angle, but that is beyond the scope of this study, and a design value 
rQ - 4 cm will henceforth be assumed here. Then, the RMS requirements for global tilt, piston, and 
subaperture tilt correction would be 3.30 prad, 42.4 rad, and 8.53 prad respectively for a 12 m aperture 
with 4 cm elements. Using the 2.5X multiplier, these would translate into P-P values of 4. Ip for edge 
stroke of a 1 m global tilt corrector (12X magnification), 13.5p for piston, and 0.05p for additional edge 
stroke of the piston to correct for subaperture tilt. 
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Spatial variations translate into temporal ones when the turbulence moves across the optical path (or vice 
versa), and the power spectrum measures the variance per unit frequency (radian 2 /Hz) of the wavefront there 
Greenwood and Fried (Ref. 2) showed that the power spectra of tilt and piston could be represented respectively by 

F a (f) = 2.079 (A/D) 2 (D/r 0 )5/3 { fc n 2 /f 0 Ha(f /f o) dz}/{ (C n 2 dz } 


F$ (f) = 2.079 (D/r 0 ) 5/3 { Jc n 2 /f 0 H 0 (f/f o ) dz } / { /C n 2 dz } 


where f 0 = V/icD is a function of the local crosswind velocity V(z), and the H’s are complicated integrals over the 
spatial region of interest Their results showed that the segment diameter d has an influence only at high frequencies 
for which f > fg D/d, whereas segment location and overall tilt both govern the low-frequency behavior for f * f 0 , 
overall tilt removal causing a significant reduction of the power densities there. With that in mind, they derived 
"single-point" asymptotic solutions for infinitesimally small segments (d = 0), both with and without prior global tilt 
removal, and then fit the data for finite element sizes with low-pass filter functions to approximately describe the 
rolloff at high frequencies due to integration over the extended transit times associated with finite d. The latter were 
found to be 


G a (y) = 1 . y < 0 332 G (p (y) = 1 , y < 0 366 

= 1.12- 0.361 y =1.26- 0.712 y 

= 0 , y > 3.10 =0 , y > 1.77 

where y = f/f 0 d/D = irfd/V is a function of the controlled segment size. In order of application for the case where 
global tilt is first removed, , the appropriate asymptotic single-point solution for "residual" tilt correction with the full 
aperture D being controlled is 


H a (x) = (2«)- 2 x- 2/ 3 , x = f/f 0 


corresponding to global tilt left in, the appropriate solution for residual piston correction of the full aperture D using a 
control size d is 


H<p(x) - 0.01037 x4* f or x < 2.216 (averaged over the region) 


= 0.25 x ‘ 8/3 


forx > 2.216 


with global tilt removed, and the solution for final residual tilt correction by the subelements would be the same as for 
the overall tilt correction but with the subaperture size d used throughout as for the Noll variances. Note that the 
high-frequency tilt cutoff is essentially the inverse of the transit time across the control element as might be expected, 
and that the rapid F 07 J fall-off of the phase spectrum will cause the subelement size to have little effect on the overall 
residual phase variance (except for the high-frequency cutoff itself) when d « D. The full-range frequency integrals 


I F a (f) df = 1 81 (D/r 0 ) 5/3 (A/nD) 2 and J F <p (f) df = 0.142 (D/r 0 ) 5/3 


of the resulting expressions are in remarkably good agreement with Noll's exact calculations and serve to validate the 
approximations that were made. 
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Because of their dependence on fp. the H’s vary along the path when the crosswind velocity does, but, for 
the sake of generality, Greenwood and Fried chose not to perform the path integrations in their paper, rather 
ascribing their results to turbulence in a single "layer" and leaving it to the user to extend the analysis as needed to 
variable crosswind scenarios (Greenwood later performed the path integrations for the phase spectra for 
representative conditions in Ref. 3.) The path integration is absolutely necessary for SELENE, especially with the 
high slew rates associated with low-orbit applications, and this is facilitated by casting the spectra equations in the 
forms 


F a (f) / {(A/D) 2 (D/r 0 ) 5/3 D/V ref } = 2.079 [ (fC n 2 /v 1/3 dz} + {1 . 12/C n 2 /v 1/3 (1 - p/v) dz} ] / {26 96 p 273 |C n 2 dz} 

for v/p > 9.46 for 1 < v/p < 9.46 


and F^ (f) / {(D/r Q ) 5/3 D/V ref ) = 2.079 p 473 [ {fC n 2 /v 773 dz} + {0.2475/p 4 /C n 2 v 573 dz) 1 / {6.67l/C n 2 dz) 

forv/p> 1.418 for v//3 < 1 .41 8 


where p = fD / V re f, v = V / V re f where V ref is some convenient reference crosswind velocity, and where we have 
again neglected the high-frequency phase rollover. These have been evaluated numerically and are shown in Table 
I and Fig. 1 for the Hufnagel-Valley profile 


C n 2 (h) = 5.94 10 

-23 (w/ 27)2 h 10 

exp(-h) + 2.7-10- 

16 exp(-h/1.5) + 

1.7 10" 14 exp(-10h) nr 273 

( h in km) 

f D / V rpf 

Uniform Wind 

Button Wind Profile 

Slew 




No Spike 

Spike 

No Spike 

Spike 

0.01 

1.661 

2.997 

2.896 

1.369 

1.395 


0.00067 

0.04333 

0.03971 

3.627 

3.290 

0.0316 

07710 

1.264 

1.228 

0.1884 

0 2417 


0.00312 

0.2011 

0.1843 

0.5814 

0.5278 

0.1 

0.3579 

0.2929 

0.3004 

0.02970 

0.05939 


0.01446 

0.9335 

0.8556 

0.1046 

0.09647 

0.3162 

0.1272 

0.00136 

0.01074 

0.00145 

0.01281 


0.06714 

0.08986 

0.09441 

0.01678 

0.02242 

1 

0 

0 

0 

2.3e-6 

0.00089 


0.07713 

0.00443 

0.00870 

0.00099 

0.00856 

3.162 

0 

0 

0 

1.6e-13 

7.4e-10 


0.00358 

0.00021 

0.00040 

4.6e-5 

0.00045 

10 

0 

0 

0 

0 

0 


0.00017 

9.5e-6 

1.9e-5 

2 1e-6 

2.1e-5 


Table I Tilt and Piston Power Spectra for Kolmogoroff Turbulence 


E.4 







'DAS* 


Rg. 1 Tit and Piston Power Spectra 
with Global Hit Removal 


without (W = 0) and with (W = 27) the "tropospheric 
spike” and for the important cases of constant crosswind 
(aka Greenwood-Fried), the Bufton wind profile 


V = 5 + 30 exp[-((h-9.4)/4.8) 2 ] mps 


with V re f = V max , and slew alone with V ref equal to the 
velocity at 10 km; the profiles are shown for convenience 
in Fig. 2. The normalization using r 0 is such that the 
results are independent of the zenith angle and it is only 
the shape, not the absolute magnitudes, of the C n 2 and 
velocity profiles that affect the results. 


Wind Speed (mp4) 



Recalling that the tilt spectra vanish for fD/V > 1, there is no contribution beyond fD/V re f = 1 for either the 
uniform or Bufton wind cases because each has a well-defined maximum velocity; the contribution to the Bufton tilt 
spectrum beyond fD/V ref = 1/7 is due solely to its shear component in the tropospheric layer. The tilt spectrum for 
slew persists to high frequencies because the velocity rises continuously with altitude, but it, too, eventually vanishes 
at about fD/V re f = 1 and 2 respectively as the high-altitude turbulence does at about 10 km without the spike and at 
about 20 km with it. The low-frequency spectra for the uniform and Bufton wind cases arises primarily from the 
strong turbulence in the ground layer, while that contribution is totally depressed with slew alone. The phase spectra 
shown for all the wind profiles persist to high frequencies where they decay as but, had the rolloff been 
considered (as will be done shortly), they would also eventually vanish. The particular value of the "reference 
velocity" clearly has a significant effect on the spectral requirements, and it is important to put that in perspective for 
the slew case: if we consider an overhead path to a circular orbit with no earth rotation, the reference velocity would 
be 


Orbital altitude, km: 

300 

1000 

3000 

10000 

30000 

Orbital velocity, mps: 

7720 

7340 

6510 

4930 

3310 

Reference velocity, mps: 

257 

73.5 

21 7 

4.93 

1.10, 
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so, with typical natural winds, slew will totally determine the high-frequency spectra for beaming to low orbits, but it 
will otherwise have little influence. Composite cases can be modeled with sufficient accuracy for most purposes by 
superimposing and adding the appropriate curves with the proper reference velocities, then scaling their strength by 
keeping in mind that the integrals must satisfy the full-range frequency integrals - which are totally independent of 
crosswind conditions. 

By definition, the integrals under the spectral curves are the contributions to the variances in that frequency 
integration interval, so, following Greenwood (Ref. 3), the residuals after correction up to a high-frequency cutoff f c 
for closed-loop servo systems with binary (sharp-edged) and a more realistic RC response are given by 

a 2 = /F(f)H(f/f c )df 


wnere H(t / r c ) = o tor r < r c , i for r > r c (Binary) 

H(f / fc) = f 2 / [f 2 + f c 2 I (RC) 


are the power rejection response functions; note that 
the low-frequency limit will generally be DC. Both the tilt 
and piston high-frequency rolloffs are important at the 
conditions of interest, and the resulting residuals are 
shown in Fig. 3 where, in all cases, d should be 
considered to be the size of the appropriate control 
element An interesting analytic solution can be 
obtained when the spectra F(f) can be represented by a 
power law of the form f n , and it can then be shown that 
the binary cutoff is related to the residual variance by 
F(rofc/V re f) rQfjVVfgf = (n - 1) o r 2 and that the RC cutoff 
frequency exceeds the binary value by a factor of [jr(n - 
1)/2] I sin[n(3 - n)/2] for a given residual; the cutoffs for 
this case are intrinsically independent of any mirror 
diameter for a given rg, and the cutoff ratio is 2.62 for n 
= 8/3 

Although there is little advantage to be gained 
by compensating global tilt to much less than a 
diffraction angle for almost any laser beaming 
application, the value on the ordinate of the global tilt 
plot will still have to be no more than about (rg/D)^ 2 ; 
this will usually be less than about 10‘ 2 for most short- 
wavelength ground-to-space situations, so the variance 
will decrease rapidly with relatively small Increases of 
cutoff frequency and can be made to essentially vanish 
with little cost penalty. Furthermore, except for beaming 
to low orbits, the transit time across the beam will 
generally be large enough that the cutoff frequency for 
global tilt control will be relatively small, typically just a 
few Hz, but caution is advised because the steepness of 
the curves is "a two-edged sword", and insufficient 
response can produce a dramatic Increase in centroid 
wander to the point that the beam becomes useless. 
For correction over the scale of a subaperture, d/rg will 
be close to unity for both piston and tilt, so the tilt 
ordinate should be about 0.01 or less QJ20 RMS at the 
edge) while the piston should be less than about 0.1 
(90% Strehl), and close examination of Fig. 3 then 
suggests that the two give about the same result. The 
resulting frequencies will be much larger, potentially 



fdAU 

Fig. 3 Residual Variances 
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several KHZ with low-orbit beaming, than those for global tilt because the relevant transit time will be that over the 
much smaller subaperture, and the warning about insufficient response bears repeating here because errors in this 
case will produce small-spatial-scale disturbances that scatter the energy over wide angles. It should be noted that, 
although the crosswind V 0 for off-zenith conditions is related to the zenith wind Vq by 

(V e / V 0 )2 = 1 - [sin 0/(1 + h/R E )] 2 . R E = 6370 KM, 

and that this will tend to compensate somewhat for the smaller rg there, the high-frequency requirements will be 
primarily determined by the minimum transit time, and that should be evaluated with the diameter d equal to the 
minimum rg (i.e. , the off-zenith value which defines the required subaperture size) and with the maximum velocity that 
will be encountered (i.e., that at zenith ) 

Turning again to SELENE with D = 12 m and d = rg 6 = 4 cm, it is useful to consider two cases, one 
with beaming only to high orbits where the velocity would be given by the Button profile (V^f = 35 mps), 
and the other for beaming to a 300 km orbit for which slew would be dominant (V^f = 257 mps). Then, 
using RC servo response with the spiked turbulence profile and with o 6 = 0.5 A/D, the normalized 

residual tilt variance would be (rtfD^iOfP/ A) 2 = 1.8 10r 5 , and the corresponding cutoff frequencies 

would be 2.9 HZ (fD/V^f - 1) and 39 HZ (fD/V^f - 1.8) respectively. Basing the subaperture 
requirements on - 0.1 and = A/IOd as before but with rgg = 8.5 cm (the zenith value), the 

normalized bit and phase variances of 0.035 and 0.35 give essentially the same requirements, 130 HZ 
for the high-orbit system (fd/V rer = 0. 15) and 1350 HZ for the low (fd/V re f= 0.21.) 

Low-orbit applications clearly impose extremely severe requirements on the AO system (and, of 
course, on the slew mechanism), so a shared system with a much smaller dedicated beam director - 
which would be effective at short ranges - might prove optimum for beaming to LEO with far fewer high- 
frequency elements. For example, a 2.4 meter telescope operating at 0.8p will produce a 2.4 m spot at 
3000 km range, and this will require only 3600 4 cm elements as compared to the 90000 elements in a 
12 m system - all of which would have to respond at high frequency! Although the global tilt control of 
the smaller system would have to respond to about 180 HZ, the required throws of all of the actuators 
would be about 4X smaller, and the larger system could be designed for the less-stressing high-orbit 
applications alone. 


1. R. J. Noll, "Zemike polynomials and atmospheric turbulence", JOSA66(3), pp 207-211 (1976) 

2. D. P. Greenwood and D. L. Fried, "Power spectra requirements for wave-front-compensative systems " 
JOSA 66(3), pp. 193-206 (1976) 

3. D. P. Greenwood, "Bandwidth specification for adaptive optics systems," JOSA 67(3), pp. 390-393 (1977) 
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DETERMINISTIC CONTROL OF HEXAGONAL SEGMENTED MIRRORS 


Dr. Glenn W. Zeiders, The Sirius Group 


Equations are presented for the control of hexagonal segmented minors based on the approach recently 
suggested by Enguehard and Hatfield (Ref. 1), examples are given to illustrate the technique, and comments are 
offered on its use. 


To review the basics of the AMP approach, consider an array of segments with mid-point edge 
differences 


dy = Uj - Uj - (tj + tj) k / 2, 


where lij is the centerline piston displacement of segment i and tj^ is the edge tilt displacement of the particular 
edge. We can then identify loops about each of the inner vertices through the mid-points and about any internal 
"holes", and it can be shown easily that the loops will be independent if each is minimized and passes only 
through segment boundaries that are internal to the surface. The sum fj of the edge differences around each 
such loop is independent of the piston displacements and depends only upon the neighboring tilts, so, since the 
tilts are the sole measureables with a Shack-Hartmann sensor and are the only knowns at the start of the 
computation, the fs may be regarded as constants. Seeking as AMP did to minimize the sum of the squares of 
all of the mid-point edge differences, subject to the constraints that each fj is constant for a given set of tilt 
measurements, the optimization problem can be expressed as 


(S/Sdj) { Z dj 2 - X f } = 0. 


thereby associating an element of a new vector X with each element of f, the elements of X (each associated with 
a different loop) being given by 


2dj = X Xj or 2d = M 2 X 


where the summation extends only over those loops actually crossing the particular edge. Substituting the latter 
in the equations for f, we obtain 


2f = M 1 X, 


so that the "optimum" edge displacements are related to the tilts by the linear transformation 

d = (M 2 M 1 - 1 ) f, 


which then directly gives the piston displacements through the first set of equations. The two M matrices - and, 
most importantly, the matrix product M 2 • - are totally predetermined by the array geometry and do not 

depend on the atmospheric distortion, so the transformation need be accomplished only once and can be realized 
in hard-wired circuitry. 
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To apply the approach to hexagonal segmented mirrors, consider the nomenclature defined by the fiaure 
at the right and by 


d i,j A = U i,j+1 - U i,j - fyj+l + tjj)A 1 2 
d i,jB = u l+1 1 J- u i,j-( t i+1,j + t i,j)B /2 
d ij c = u i+1,j-1 ' u i,j ' Ctj+i j_i + tj j) c 12 

where the tilts in the triangular frame are related to Cartesian ones by 
t A = <t x + V3ty)/2, % = t x , tc = ( t x - V3 ty ) / 2 



Thus, traversing the triangular loops In the clockwise direction, the constraint relations become 


f i,j d i,j A - d i,J B + d i,J+1 c = - (‘ij+l + t i,j)A 1 2 + (*1+1, j + *i,j)B 1 2 * (*i+1 J + ^.j+lJc 1 2 
= I ^ (t|+i j - tjpy + (tj+ij + tjj - 2 tjj +1 ) x ] /4 

f i,j = - d l+1,i-1 A + d i,j b - d l,j c * (*i+1J + t i+l 1 j-l)A /2 -(*i+lj + t i,j)B /2 + ( t i+1,j-l + t i,j)c /2 
= I * (*i+1,j - tjj)y - (tj+1 j +t, j - 2 1( +1 j. l)x ] / 4, 


and the result of the differentiation is 


2 d i j A = ^ij'^’i-lj+l 
2 d i,j B * - j + *'ij 

2d i,jc * Hj-i'Hi 


considering only triangular loops about internal vertices; X terms that correspond to “open" loops at the external 
boundaries will vanish, and terms at inner boundaries will be replaced by ones for the longer loops used to 
encircle holes in the array. Substituting the above in the equations for f and f , we obtain 


2 *i,j ~ 3 \j" ^'i-lj+1 ' *'i,j+1 ’ *‘ij 
2f i,j = 3 *'ij * *ij-1 • *i,j - *i+1,j-1> 

each of which effectively averages the primary triangular path with its three nearest neighbors - and enables 
relatively easy construction of the and 1*2 matrices as shown below. If segment i,j were to be removed, the 



U ^i+1 i-l 1 i i+1 ^il+1 1 

2f|i 

3 -V !!!!! ' -V ' Vi" . VV 

-1 -1 -1 3 

2diiA 
2di j r 
2dj i 

c 

~~ ~ ~ ~ ~ i rr: ^ ~ 7 ~ ~ 

-1 1 

-1 1 . 




equations for fj_i j. H-ij+i Jjj. fjj, , and fj_-j j would be replaced with the single equation 

2 f H = 6 A. H - *V2,j+1 ' *i-1 j+1 - *‘i,j+1 ' *i+1,j - ^’ij-1 - *i-1 J-1. 


for the six-sided loop around the hole, the edge differences dj j A . dj i B , dj i c , d t _ 1 j A , d ( . 1 j B , and d|.-j j c would 
vanish, and those for dj.-| j A , dj_i j + .j B , dj j +1 c , dj +1 j_-j A , djj_’-j B , and dj.^’ j c would be replaced with 


2 d M j A = ‘ 2 dj-ij+1 B = " *i-1,j+1 2 dj j + -j c = * ^'ij+1 

2 d i+1,j-1 A 3 X i+1,j-1 * X H 2 d i,j-1 B = X ’i,j-1 * *H 2 d M,j C 3 *i-1 j-1 ' *H 


A representative set of loops and the corresponding Mi and M 2 matrices are shown in Fig. 1 for a 19- 
element filled array, and the solution matnces M^and M 2 - M^** are given in Fig. la. Only five types of mating 
edges are possible with a 19-element filled array, and these are shown in Fig. la with the corresponding 
weighting functions from M 2 • M^’ 1 ; as an aid to visualizing the extent of the interactions, loops whose weights 
exceed 10% of the maximum are shaded, and those which exceed 20% are the darker. Equivalent results for the 
same array with a missing central element are shown in Figs. 2 and 2a for comparison (there are only four types 
of mating edges with that geometry.) Although the input matrices are quite sparse, the solution ones are not, and 
the results show that the influence generally extends well beyond the nearest neighbors, suggesting that 
convergence problems might be encountered with an iterative solution. Note in particular that the missing 
segment actually causes the influence to be more widely spread, causing further problems for iteration while 
having no negative impact on the deterministic solution. 

This advantage of the AMP approach is that it is deterministic and not iterative, but its success depends 
upon the ability to establish a reference surface from which piston displacements can be measured. Although the 
predicted edge differences can in theory be used to individually adjust the segments (or each be combined with 
the local average tilt to give the piston differential across the gap), the differences are likely to be so small that 
measurement errors will accumulate excessively, and they should be used instead in groups to specify piston 
displacements from a common reference piston - not each with reference to a nearest neighbor. The flexible 
telescope structure being pursed for SELENE appears to preclude an overall mechanical reference for piston 
displacement, and the atmospheric distortions from the "beacon" prevent it from being an optical reference; 
similarly, the edge measurement techniques that have been proposed may not be sufficiently accurate to prevent 
an unacceptable error buildup across the beam if edge displacements were to be used throughout. However, a 
"hybrid" system (woofer/ midrange/tweeter) using the flexible substructure with rigid multi-element segments, 
either interferometric or edge-referenced, might provide a workable compromise. Otherwise it appears that an 
"invention'" will be required to provide the required reference. 

A second key question of the importance of the analysis is, of course, whether minimization of the mid- 
point edge differences really represents the optimum condition. It dearty does not offer the precision of a 
wavefront-matching approach such as that claimed by JPL, but ifs a physically realistic one, and it deserves to be 
tested, probably via a Zemike analysis against a more rigorous approach, to see which offers the greatest cost 
benefit It should be kept in mind that AMP made no claim that the mid-point matching condition was indeed the 
best one, but they chose to use it since it should be near-optimum and since it had already been accepted as an 
optimizing condition by KAMAN and others. Their approach could probably be modified to accept a better 
premise should one be offered. 


1. S. Enguehard and B. Hatfield, "Minimizing Edge Misalignments for Segmented Mirror Adaptive Optics," 
AMP-92-10 (1992) 
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Figure 1 19-Element Hexagonal Array 
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1 9-Element Hexagonal Array Solution 
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Figure 2 19-Element Hexagonal Array 
with Missing Central Segment 
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Figure 2a Solution for 19-Element Hexagonal 
Array with Missing Central Segment 
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SOME ASPECTS OF PHOTOVOLTAIC POWER 
CONVERSION WITH LASER ILLUMINATION 


Dr. Glenn W. Zeiders, The Siriue Group 


The current-voltage characteristic of a p-n junction photovoltaic converter is given by 


l='«- UeWVi) 


where 1^ is the short-circuit current when the cell is illuminated, and l 0 and V 0 = nkT/e are characteristics 
of the cell. The parameter n is ideally unity but can be as much as three for real systems, kT/e = 26 mV 
@ room temperature, and the open-circuit voltage = V 0 /n(1 + l^/lp) is about 550 mV for Si at one sun 
and rises linearly with bandgap, so l^/l 0 is generally quite large for cases of interest. Since photon flux 
is effectively converted to electron flow, the short-circuit current l K is proportional to the incident radiant 
intensity. The simplest realistic equivalent circuit for a cell is a current source with a series resistance 
R,, so the overall circuit behavior for an array of identical cells with uniform illumination is described by 


X=ln{l K /l 0 + 1-Y}-Y (I.R./VJ 


where X = V/NV 0 and Y = l/l 0 for series-connection 
of N cells, X = V/V 0 and Y = l/NI # for parallel- 
connection, and where V and I are the system 
output voltage and current The conversion efficiency 
is proportional to X Y = IV / V 0 L,N, and is realized 
for 

d(XY)/dY = ln{l lc /l 0 +1 - Y}-Y/{ l^/l + 1 -Y} 
-2Y(l o R,/V o ) = 0, 



Rg.2 Photovoltaic Response 
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the result of which is plotted in Fig. 1 . The optimized 
"efficiency" rises slowly at first because of the 
logarithmic variation of V^, but then falls abruptly 
due to l 2 R § ohmic heating in the series resistance. 
Because the optimum current is proportional to the 
l, c and the optimum voltage is proportional to V^, 
foe optimum load resistance R = V / 1 varies roughly 
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inversely with l lc , so the series resistance becomes more critical at high illumination levels. The effect of 
the total resistance (load and internal) is shown in Fig. 2, from which it can be seen that a very serious 
performance penalty is paid if the resistance is not within at least a factor of two of the optimum value. 
This can easily be handled with solar flux alone, but the required accuracy and dynamic range may be 
quite difficult to achieve if the system must also accommodate large peak currents from laser pulses. 

The problem is further complicated with laser illumination because, unless the array is flood- 
loaded with the beam, the illumination will vary both spatially (due to the beam intensity profile) and 
temporally (due to beam wander from system jitter and atmospheric turbulence along the path.) It has 
been suggested that cells be connected in series to provide high voltage at reduced load current, thereby 
decreasing the series-resistance loss, but that does not work well when the irradiation varies appreciably 
over that portion of the array. The impact of the type of connection is shown quite effectively in Fig. 3 
with a simple two-cell example for which series and 
parallel connections are compared with both 
matched and unmatched cells. Internal resistances 
are neglected in the figure, and it is assumed in the 
latter case that the weaker cell has half the 
illumination of the other (l 8C /l 0 =210® for the higher.) 

Load lines are shown for optimum operation with the 
matched cells, and, while the total incident power 
decreases to 75% for the unmatched case, the 
electrical power with series connection drops to 31% 
whereas it falls only to 71% with parallel connection. 

The reason for the nearly constant efficiency with 
parallel connection is that the open-circuit voltage will 
be changed little by mismatching (again because of 
the weak logarithmic dependence on El sc ), and the 
total current will continue to be near the sum of the 
individual short-circuit currents. On the other hand, 
the maximum current with series connection will be 
no more than about the short-circuit current of the v/v 0 

weakest cell, and that cell's open circuit voltage will „ o _ 

depend on it's short-circuit current alone, so the ” 0 * “ Ltt©Ct Of C©ll Connoctlon 

efficiency can be dramatically reduced by mismatch- 
ing as found in the example. In the worst case of complete "shadowing", weakly-irradiated cells in 
series not only decrease performance, but they can also act as dead shorts capable of causing severe 
damage through "hot-spot" formation. Diodes can be used to protect against such damage, but they do 
not overcome the performance loss, and they add weight and complexity to the system; some protection 
is afforded as well by leakage currents arising from defects and edge effects in the cells, but their effect 
is generally unimportant otherwise. The net effect of this is that series connection can — and probably 
should -- be used over regions small compared to the instantaneous size of the main beam in order to 
obtain the desired voltage while reducing series resistance loss, but these elements should then be 
either independent or connected in parallel to maximize efficiency over regions with significant intensity 
variation. 

The response to a directed beam (especially one from the free-electron laser of interest for 
power beaming) can be modeled by assuming a gaussian distribution of intensity and, in view of the 
above results, a system of individual elements each with its own load resistance. Neglecting internal 
losses, the converted power distribution is 


P’,r, = 2 * JJ 2r *he„ 

Case I: R = constant or Case II: l/l 0 = (V^R) fn (V^R) for local optimization, 
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and where R is the resistance per unit area of 
the array. The total power output for the two 
cases (again normalized with respect to l K to 
represent "efficiency") is shown in Fig. 4; the 
maximum power output with uniform resistive 
loading and gaussian illumination is achieved if 
the load is optimized for about 55% of the peak 
short-circuit current, and the performance 
penalty for not using locally-optimized loading is 
27% at t sc /l 0 =10® and improves to 15% at 10 3 . 
This rather nominal performance penalty with 
fixed loading may, indeed, prove to be a very 
small price to pay when one considers that the 
locally-optimized loading would have to be 
"active" at relatively high frequencies to account 
for centroid motion; the constant-resistance 



Rg>5 Distribution of Optimized 
Response to Gaussian Irradiation 



"passive" load might have to be variable as well 
(particularly for space vehicles) to account for changes 
of peak intensity, but those variations would be applied 
on a global scale at much lower frequencies. Fig 5 
shows the distribution of output power density and 
included power for (l^/lo)^ = 10 6 for the optimum 
fixed resistance as well as for locally-optimized 
("active") loading. The value of the fixed resistance is 
higher than optimum in the high-intensity central 
region and lower than optimum in the wings, but the 
penalty is rather nominal as noted earlier, and the 
output distribution is somewhat compressed. For 
example, 86.5% of the incident power is included 
within the 1/e 2 radius r = 2a, whereas, for the case 
shown in Fig. 5, that same fraction of the output power 
is included within r/a =1.89 with local optimization and 
only 1.55 with the optimum constant load, thereby 
providing the possibility of a smaller receiver. 
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